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We analize the stability of a class of thin-shell wormholes with spherical symmetry 
evolving in flat FRW spacetimes. The wormholes considered here are supported at 
■ the throat by a perfect fluid with equation of state V = wa and have a physical 

, radius equal to aR, where a is a time-dependent function describing the dynamics 
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is done by means of the stability analysis of dynamic systems. 
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Introduction 



Wormhole physics has become a popular research topic since the classical paper by Morris 
and Thorne.^ Early work was reviewed in the book of Visser ^ and there is an extensive 
recent review by Lobo.'^ The class of thin-shell wormholes provides a wide collection of 
examples for the study of traversable Lorentzian wormholes. Visser ^'^ studied wormholes 
I in Schwarzschild and Reissner-Nordstrom backgrounds using for the first time the cut-and- 
past technique. The stability of static and dynamic wormholes has been examined either by 
choosing specific equations of state or by considering a linearized stability analysis around 
a static solution. The literature is extremely extensive but for some examples see Refs. 
6-19 and references therein. In this work we analyze the stability of a class of thin-shell 
wormholes with spherical symmetry evolving in a flat Friedmann-Robertson- Walker (FRW) 
cosmological background. Dynamic Lorentzian wormholes connecting FRW spacetimes were 
first studied by Hochberg and Kephart ^ who discussed a possible resolution of the horizon 
problem using a network of evolving wormholes present in the early universe. We utilize the 



* Electronic address: |lacamera@ge.infn.it| 



2 



wormhole equations of motion given in Ref. 7 and assume that the matter at the shell is a 
perfect fluid with equation of state V — wa. We shall obtain an evolution equation for the 
physical radius aR of the wormhole, where a is a time-dependent function describing the 

dynamics of the throat and R is the background scale factor. The study of stability is done 
considering the wormhole physical radii as dynamic systems and analyzing the behavior of 
their steady states with respect to a set of control parameters affecting the evolution. 



2. Thin-shell wormholes in cosmology 

The dynamic wormholes we shall consider here result from surgically grafting two FRW 
spacetimes as proposed by Hochberg and Kephart ^ following the method used by Visser 
^'^ for constructing wormholes by cutting and pasting two manifolds to form a geodesically 
complete new manifold with a shell placed in the joining surface. The construction of Ref. 
7, which we now briefly recall, starts taking two copies of the FRW solution 

\1 — kr'^ ) 

and then removes from every one of them an identical four-dimensional region of the form 
E± — {r± < a} where a = a{t) is a time- dependent radius. The resulting manifold contains 
two disjoint boundaries E-t = {r± — a} and by the identification = E_ = E one obtains 
two FRW spacetimes connected by a traversable wormhole with spherical symmetry whose 
throat is located on their mutual boundary E = {/(r, t) — r — a{t) — 0}. The corresponding 
geometry can be analyzed using the Israel-Lanczos-Sen 20-22 "thin-shell" formalism of Gen- 
eral Relativity. The wormhole throat is a timelike hypersurface with interior coordinates 
^* = {i9, ^p,t), t being the throat proper time. The position of the throat in the background 
FRW spacetime is given by — {a{t),'d,if,t). The line element (1) with the curvature 
term k neglected is well justified to describe not only, as in Ref. 7, the early stages of 
expansion but also the universe today on large scales, so we shall put A; = in it. Moreover 
another realization of an evolving wormhole was obtained by imagining its asymptotically 
parts as constituting flat FRW spacetimes. ^ The tangent vectors to S have the components 
^('j-) |± = dX^/d^'^ and the induced metric gij = e'^--^ e'^j^ g^y \± on the junction hypersurface 
gives the line element 

ds^ = a^i?^ ((i^2 ^ g.^2 ^^^2-) _ ^^2 (2) 
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having set 



dt 



+ 



dT 



a{t)] R^it) 



—1 the coefficient of dr"^. Notice that the 



dt 



physical radius of the wormhole is aR. Using a' = a—, where the prime and the dot 

dr 

represent the derivatives with respect to r and t respectively, one obtains 



dt 



1 



(3) 



dT - {dRf 

therefore it will be possible in the following to eliminate one of the two time parameters r 
and t in favour of the other. The second fundamental form (extrinsic curvature) is given by 



•9^^^ ^-v± dX^ 



a/3 



(4) 



where is the unit normal to E: 



± 



Of df 



dX" dXP 



-1/2 



df 
dX^ 



(5) 



With the definitions Kij — — K^^j the Lanczos equations follow from the Einstein field 
equations for the hypersurface, and are given by 



(6) 



where Sj is the surface stress-energy tensor on E. In the case of spherical symmetry consid- 
erable simplifications occur, namely K^j — diag («;|, k^, k^) , and the surface stress-energy 
tensor may be written in terms of the surface pressure V and the surface energy density a 
as Sj = diag {V, V, — a). The Lanczos equations then reduce to 



a — — — K,< 



(7) 
(8) 



where 



(iR\ 



aR 



a"R 



+ 2a'R 



(9) 
(10) 



Using Eq. (3), it is possible to substitute derivatives with respect to the proper time on the 
shell r with derivatives with respect to the cosmic time t, so the components of the tensor 
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Sj are writtem as 



,21-1 „^ , ^. A , 1 



[1 - - (ciR) + 2dR + ^ 

7? ^ ^ <U (jJi (II) 

47r - (ai?)2 

di?+ — 

a = , (12) 

27r VI- 

and the conservation identity = 0, which is not independent but can be obtained from 
the field equations, provides the simple relationship 

If one chooses a particular equation of state, in the form V — V{cr), then the conservation 
equation can be formally integrated obtaining 

da 



The choice V — V{a), because of Eqs. (11) and (12), will set a bond between a and it!. In 
this paper we shall consider the case 

V = w(T (15) 

where w is an arbitrary real constant. Therefore the conservation equation (14) and the 
equation of state (15) become respectively 

- TD \ 2(i+«') 
0,0 Ko 



[1- imf]-' |(<ifl) + 2(l + io)afl+^iip^=0 (17) 

where the subscript o refers to quantities calculated at some initial time t — Iq. We have 
four unknowns: a, i?, V and cr, but only three independent equations: (11), (12) and (15). 
The standard stability method based on the definition of a potential ^ is not applicable here 
because from those equations it does not appear possible to obtain explicitly the derivative 
with respect to the time of the physical wormhole radius aR as a function of aR. Thus, to 
close the system, we need to have one more equation. A suitable fourth equation must be 
chosen in such a way to facilitate the stability analysis for the examined class of wormholes. 
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3. Stability analysis 

In this section, to analyze the stabihty of the wormhole physical radius we substitute again, 
using Eq. (3), derivatives of aR with respect to r with derivatives with respect to t and it 
may be checked that at equilibrium, where —{aR) and —{aR) vanish, the final result with 
regard to the stability does not change. We consider the wormhole radius aR as a dynamical 
system (see, e.g., Ref. 23) with an evolution equation of the form 

{aR) = F{aR, A) (18) 

where F is a differentiable function acting on the spacetime where aR is defined and A 
denotes a set of control parameters affecting the evolution. The value of aR for which 
F{aR, A) = is a steady state point. The stabihty of the equihbrium will be studied 
considering small perturbations away from the above fixed point and determining, by means 
of a linear analysis, the tendency of the perturbations to grow or to decay in time. Now 

^ {aR) ^aR + ait (19) 

(Jjv 

so we must know both aR and aR as a function of aR. Here it is worth noticing that when 
w — —1/2 the equation of state (17) becomes 

[1 - {dR)'^] ^{ciR) + dR = (20) 
at 

which is satisfied either when the throat function a is constant, and so d = 0, or when a is 
time-dependent, and d is given integrating Eq. (20): 




(doi?o) 
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-1/2 



(21) 



(aoi^o)' 

When d = 0, the wormhole radius aR simply varies in direct proportion to the scale factor 
R. When d 7^ 0, which is the case treated in this paper, the behavior of aR can be studied 
if the scale factor R is known. When w ^ —1/2, to choose a suitable fourth equation we 
consider a specific class of simple wormhole solutions corresponding to the choice 

[l - {aRf\ ^{aR) -2waR = Q (22) 

which generalizes Eq. (20) to the other values of w, and the equation of state (17) becomes 

{l + 2w) (2a^+^) =0 . (23) 
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In the following we shall consider values of w different from the value w — — 1/2 which was 
briefly discussed above, so Eq. (23) gives 

2 di? + ^ = . (24) 
an 

Then we shall treat a particular class of wormholes for which a and R must have opposite 
signs and therefore when one is increasing the other is decreasing. Now we substitute 
aR = — l/(2di?) into Eq. (19) which becomes 

|(afi) = dfi-^ (25) 

therefore the equilibrium is reached when (aR)^ = 1/2. Moreover the components of the 
tensor Sj now become 

w 

V ^ , (26) 

47rai?^l - (di?)2 

a ^ / (27) 

47iaR^/l - (di?)2 

so the weak energy condition is violated, while the null energy condition is satisfied when 
w < —1. Finally, to obtain aR as a function of aR we equate Eqs. (16) and (27) for the 
surface density a and have 

l + 47raoi?o<7oVl- («^)' ^ =0 . (28) 

Evaluating the previous equation at the time t = to it results 



1 + 47raoi?o(7o Vl - (ooi^o)' = . (29) 
Then Eq. (28) takes the form 



{dRf [ aR 



1 - (do-Ro)^ \aoRQ 



l+2w 



(30) 



and therefore 



aR ^ 



|di?| = ^/l- [l-(doW](^^J . (31) 

Then, when w — Q both di?, by Eq. (22), and aR are constant. Now we can write Eq. (25) 
in the form of Eq. (18): 

4- (aR) = sign[d] ° (32) 
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where sign[d] = a/\a\. The right-hand side vanishes when 

a.R. = aoi^o [2 [1 - [a^Rof]] -"^^^'^^^^^ (33) 

so a^,R^, is the static solution which gives the radius of the throat at the equihbrium. If 
[aQRof — 1/2, then the static solution is a^RQ and does not depend on w. The role of 
perturbation of stability will be expressed by setting in Eq. (32) 

aR = a^R^ + i{t) (34) 

where |,^| <^ a-t^R-t represents the disturbance. We are interested in linear analysis so higher 
order terms will be neglected in the Taylor expansion of the evolution equation (32) which 
yields at first order 

^ e = - sign[d] V2 (1 + 2w) (35) 



Therefore we can conclude, at least until the nonlinear region is reached, that : 
(i) if (l-l-2'u;) sign [d] > 0, then the perturbation ^(t) decays exponentially and the equilibrium 
is stable; (ii) if (1 + 2w) sign [d] < 0, then the perturbation ^{t) grows exponentially and the 
equilibrium is unstable. Here d and R have opposite signs, so in an expanding universe the 
equilibrium is stable when w < —1/2. Finally we recall that, having chosen the equation of 
state in the form V = V{a), the functions a and R are related. This can be seen explicitly 
integrating first Eq. (22): 

and then equating, in the case w ^ —1/2, Eqs. (31) and (36). The result is 



O RO I / . „ s On / . „ s O / -R ^ 



-l/(2(l+2w) 

[l-(doRor] + {aoRor[^] } (37) 



4. Conclusion 

In this paper we have considered a particular class of thin-shell wormholes evolving in a 
flat FRW cosmological background and analized their stability to linearized perturbations 
around static solutions using the stability analysis of dynamical systems. The physical 
radius of the wormholes is equal to aR, where a is a time-dependent function describing 
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the dynamics of the throat and R is the background scale factor. We have considered the 
equation of state V = wa and found that the factor w is crucial for determining the stability 
of the equilibrium. Wormholes supported by other equations of state than the one used in 
this paper can be used to generate other families of additional solutions. 
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